Radiative and Collisional Jet Energy Loss in the Quark-Gluon Plasma at RHIC 
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We calculate and compare bremsstrahlung and collisional energy loss of hard partons traversing 
a quark-gluon plasma. Our treatment of both processes is complete at leading order in the coupling 
and accounts for the probabilistic nature of the jet energy loss. We find that the nuclear modification 
factor Raa for neutral tv° production in heavy ion collisions is sensitive to the inclusion of collisional 
and radiative energy loss contributions while the averaged energy loss only slightly increases if 
collisional energy loss is included for parent parton energies E S> T. These results are important 
for the understanding of jet quenching in Au+Au collisions at 200 AGeV at RHIC. Comparison 
with data is performed applying the energy loss calculation to a relativistic ideal (3+l)-dimensional 
hydrodynamic description of the thermalized medium formed at RHIC. 



Introduction - Relativistic heavy ion collisions are de- 
signed to produce and study strongly interacting mat- 
ter at high temperatures and densities. Experiments at 
the Relativistic Heavy Ion Collider (RHIC) have demon- 
strated that high pt hadrons in central A + A collisions 
are significantly suppressed in comparison with those in 
binaryn + p collisions, scaled to nucleus-nucleus colli- 
sions [3, 0, H| . This result has been referred to as jet 
quenching and has been attributed to the energy loss of 
hard pr partons due to induced gluon bremsstrahlung 
in a hot quark-gluon plasma. Bremsstrahlung energy 
loss has been calculated in several theoretical formalisms 
before 0, d, 0, 0, 0, 0] ■ Recently such bremsstrahlung 
calculations were implemented in models employing rel- 
ativistic ideal (3+l)-dimensional hydrodynamics in or- 
der to calculate the nuclear modification factor Raa of 
neutral pions at RHIC [l^, El, [12] ■ Early estimates of 
the collisional energy loss which used asymptotic argu- 
ments indicated that the radiative energy loss is much 
larger than elastic energy loss Zakarhov compared 
radiative energy loss in the light-cone path integral ap- 
proach and collisional energy loss employing the Bjorken 
method and concluded collisional energy loss is relatively 
small in comparison to the radiative one Q • Renk derives 
phenomenological limits on radiative vs. collisional en- 
ergy loss by considering quadratic vs. linear pathlength 
dependence and concludes that any elastic energy loss 
component has to be small [l4 ]. In contrast, Mustafa 
and Thoma find that collisional energy loss has a signif- 
icant influence on jet quenching [TBI , [la ] . Recent studies 
by Gyulassy and collaborators also point in this direc- 
tion, sec e.g. fi~7l.ll8l]. 

The purpose of this study is to consistently incorporate 
collisional and radiative energy loss in the same formal- 
ism and to employ this formalism in a realistic descrip- 
tion of energy loss of hard px leading partons in the soft 
nuclear medium as described by (3+l)-dimcnsional hy- 
drodynamics in 200 AGeV Au+Au collisions at RHIC. 

We will emphasize three points (the first two of which 
have been elucidated earlier [l9|| for bremsstrahlung en- 
ergy loss). First, in many previous approaches the aver- 
age energy loss is computed and applied to the primary 



partons. Bremsstrahlung and collisional energy loss are 
not well described by a (path length dependent) average 
energy loss alone. Bremsstrahlung energy loss is domi- 
nated by hard emissions. Therefore, if a sample of par- 
tons initially has the same energy, then after traversing 
some pathlength of the medium, the distribution of final 
energies will be in general broad and not sharply cen- 
tered at an average energy. This will be illustrated in 
Fig. [TJ While the average energy loss has some value 
in judging the importance for observable consequences 
in jet-quenching, the evolution of the probability den- 
sity distributions of partons until fragmentation is the 
decisive quantity for such studies. To account for this 
we directly evolve the spectrum of partons as they un- 
dergo bremsstrahlung and collisional energy loss. Sec- 
ond, radiative energy loss depends on a coherence ef- 
fect: the Landau-Pomeranchuk-Migdal (LPM) suppres- 
sion. While some approaches take the LPM effect as a 
parametrically large suppression, this is only true when 
the parent parton and the emitted gluon arc highly en- 
ergetic, E'parton, -Egiuon ^ T. For small radiated gluon 
energies E g \ ucm <C -Eparton the LPM suppression can be 
significantly less, and those bremsstrahlung events are of 
significant importance in jet quenching [10] due to the 
steeply falling initial parent parton spectrum. We there- 
fore employ the Arnold, Moore and Yaffe (AMY) formal- 
ism [§] for radiative energy loss to treat the LPM effect 
at all energies £gi U on > g s T correctly up to 0(g s ). Third, 
while there has been considerable theoretical effort to im- 
prove our understanding of jet modification in the quark- 
gluon medium, early jet quenching calculations often re- 
lied on an elementary description of the soft medium. 
Until recently most jet quenching calculations used sim- 
ple medium models only loosely constrained by the value 
of bulk observables. Previously we presented a calcula- 
tion of Raa in central and non-central collisions using 
the AMY formalism and a (3+1) dimensional hydrody- 
namical model constrained by soft observables at RHIC 
[l2| ■ Here we incorporate collisional energy loss into this 
framework. 

Brief review of the formalism - For details of the cal- 
culation of the initial distributions of hard partons in the 



2 



early stage of the collision and the subsequent propaga- 
tion through the hot and dense quark gluon plasma as 
well as the fragmentation we refer the reader to [UJ and 
references therein. We concentrate here on the incorpo- 
ration of collisional energy loss in the formalism. 

The jet momentum distribution P(E, t) = dN ^^ 
evolves in the medium according to a set of coupled 
Fokker-Planck type rate equations, which have the fol- 
lowing generic form [l2l ]: 



dP(E) 
dt 



P(E- 



^dT(E+uj,uj) 
du 



P{E) 



dT(E, uj) 



(1) 



where dT(E,oj)/duj is the transition rate for processes 
where partons of energy E lose energy uj. The u < part 
of the integration incorporates processes which increase 
a particle's energy. The radiative part of the transition 
rate is discussed in [l^, [H, H(j ■ 

Now wc must add the contribution from collisional en- 
ergy loss. Compared to radiative loss, collisional losses 
are more dominated by small energy transfers; the con- 
tribution to the mean energy loss rate dE/dt from elas- 
tic collisions is logarithmically sensitive to large energy 
transfers, while the radiative contribution is power-law 
dominated by large radiations. Therefore it should be 
an adequate procedure to approximate clastic collisions 
by a mean energy loss, provided we include a momentum 
diffusion term as dictated by detailed balance. 

The leading order mean collisional energy loss rate is 
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where / is the thermal distribution of the medium par- 
tons. I M 1 2 is the ^-channel scattering matrix element 
squared calculated in leading order, and g^ is the de- 
generacy factor for the initial thermal partons. Eq. @ 
is infrared logarithmic divergent, screened by plasma ef- 
fects which are incorporated by including hard thermal 
loop corrections for soft momenta ~ gT. The resulting 
differential energy loss dE/dt\ a b for the scattering of a 
light hard parton a off a soft parton b are [28[ : 
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where c b = - lE + C'( 2 )/C(2), c f = c b + ln(2), and c s « 

27ra s T 2 (l + n f /6) is 



T. 66246 are constants and to 2 



the thermal gluon mass [21(. 



We can incorporate these dE/dt results in Eq. (JTJ 
by introducing the drag term, (dE/dt)dP(E)/dE, and 
the diffusion term, T(dE/dt)d 2 P(E)/dE 2 . We discretize 
Eq. ||T]), such that J du> — > Aui ,_„ a , . , and 



T{E + Aw, Aw) 
T(E,-Aw) 



(l + / B (Aw))(Aw)- 1 d^M, 
/b(Aw)(Aw)" 1 dE/dt, (4) 



which yields the right energy loss rate and preserves de- 
tailed balance. 

Results - In order to illustrate how collisional and ra- 
diative energy loss influence the time evolution of the 
leading parton distributions, we first consider a static in- 
finite medium with T = 400 MeV and a s = 1/3 and 
an initial single light quark jet of energy Ei = 16 GeV 
propagating through it. 
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FIG. 1: (Color online) The evolution of a quark jet with initial 
energy Ei = 16 GeV propagating through a static medium of 
temperature T — 400 MeV, where the vertical lines repre- 
sent the values of mean energy related to the corresponding 
distributions. 

In Fig. [1] we compare the evolution of the jet par- 
ton distribution P(E, t) under three different approxima- 
tions: (1) with only collisional energy loss, (2) with only 
radiative energy loss (already calculated in [19| ) , and (3) 
with both energy loss mechanisms. The first moment 
in energy of these distributions defines the mean energy 
(indicated by vertical lines) and indicates the average en- 
ergy loss. The figure indicates as expected that radiation 
lead to a larger mean energy loss than with elastic col- 
lisions only. As pointed out earlier, small differences in 
the average energy loss do not necessarily imply small 
differences in the parton distributions. While the time 
evolution of P(E, t) in case (3) resembles qualitatively 
the case (2) in which only radiative energy loss has been 
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considered, quantitative differences especiaily at energies 
closer to Ei can be large. 

We will model the behavior of the quark-gluon medium 
using relativistic fluid dynamics, which has been shown 
to give a good description of bulk properties at RHIC. 
We use a fully (3+1) dimensional hydrodynamical model 
for the description of 200 AGeV Au+Au collisions at 
RHIC [22j |. The initial momentum distribution of jets 
dN^ /d 2 p 3 T dy\i n i is computed from pQCD in the factor- 
ization formalism, for details see r23l. I24T] . The final 
hadron spectrum dN h /d 2 p T dy at high Pt is obtained by 
the fragmentation of jets in the vacuum after their pass- 
ing through the (3+1) dimensional expanding medium 



dN h 
d 2 p T dy 



J2 J d 2 r ± V(b,r ± ) J ^lD h/j (z,,Q F ) 



dN j (b,f±) 



d 2 pLdy 



(5) 



Bn 



where dN 3 '(b,f±) / ' d 2 p' T dy\f ln is the final momentum dis- 
tribution of the jet initially created at transverse position 
f± after passing through the medium. This distribution 
is calculated for every specific path through the medium 
by solving Eq.{T]) incorporating collisional and radiative 
energy loss. The fragmentation function D h /j(zj 7 Qp) 
[25l | gives the average multiplicity of the hadron h with 
momentum fraction Zj = px/p ] T produced from a jet j at 
fragmentation scale Qp. V(b, r±) is the initial jet den- 
sity distribution at the transverse position f± in colli- 
sions with impact parameter b. For further details see 
[l2| where radiative energy loss has been studied in an 
analogous manner. 

The final hadron spectrum directly enters the calcu- 
lation of the nuclear modification factor Raa which is 
defined as the ratio of the hadron yield in A+A colli- 
sions to that in p+p interactions scaled by the number 
of binary collisions N co \\: 



R h AA {b,p T ,y) 



1 dN h (b)/d 2 p T dy 



(6) 



Once temperature evolution is fixed by the initial con- 
ditions and subsequent evolution of the (3+1) dimen- 
sional expansion, the strong coupling constant a s is the 
only quantity which is not uniquely determined by the 
model. In this work we take its value to be constant at 
a s = 0.27, which reproduces the most central data [Hj]. 

In Fig. [2] we show the mean energy loss of quark jets 
passing through the nuclear medium created in central 
collisions (b = 2.4 fm) at RHIC, as a function of their 
initial energy Ei . In this figure, the jets are assumed 
to be created at the center and propagating along the 
in-plane direction. In agreement with [31 we find that 
the average energy loss is not strongly changed by ac- 
counting for clastic collisions. In Fig. [3] we present the 
calculation of Raa f° r neutral pions measured at mid- 
rapidity for two different impact parameters, 2.4 fm and 
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FIG. 2: (Color online) The mean energy loss of a quark jet 
with initial energy Ei passing through the nuclear medium 
created in central collisions (b = 2.4 fm) at RHIC. The jet 
starts from the center of the medium and propagates in plane. 
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FIG. 3: (Color online) Nuclear suppression factor Raa for 
neutral pions in central and mid-central collisions. Here the 
dashed curves includes only radiative energy loss, the dash- 
dotted curve includes only collisional energy loss and the solid 
curve includes both radiative and collisional energy loss. 



7.5 fm, compared with PHENIX data for the most cen- 
tral (0-5%) and mid-central (20-30%) collisions 0. We 
present Raa for purely collisional (1) and purely radia- 
tive (2) energy loss, as well as the combined case (3). 
One finds that while the shape of Raa for case (3) is not 
strongly different from case (2) that has only radiative 
energy loss, the overall magnitude of Raa changes sig- 
nificantly. We checked (comparison not shown) that the 
stronger influence on Raa stems from the differences in 
the evolution of the parton distributions in case (2) and 
(3). This has already been discussed in the static case 
(compare Fig. [IJ. The magnitude of Raa is therefore 
sensitive to the actual form of the parton distribution 
functions at fragmentation and not only to the average 
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energy loss of single partons (compare Fig. [2]). In [l2j], 
the observational implications on Raa measurements due 
to only radiative energy loss were studied. Recalculating 
Raa versus reaction plane including elastic collisions we 
found only small differences (after adjusting the coupling 
strength from a s = 0.33 to a s = 0.27) in the shape of 
Raa as a function of pt and the azimuth. 

Conclusions - We calculated collisional and radiative 
jet energy loss of hard partons in the hot and dense 
medium created at RHIC. We treated the LPM effect us- 
ing the AMY formalism Q and treated collisions using a 
drag plus diffusion term reproducing dE/dt and detailed 
balance. While we find in accordance with 0] that the 
additionally induced average energy loss due to elastic 
collisions is small in comparison to the radiative one, the 
time evolution of the parton distributions P(E, t) in both 
cases differ significantly. This is especially true for ener- 
gies close to the initial parton energy. Since the initial 



parton spectrum is steeply falling, stronger differences in 
Raa can result. We find that the inclusion of collisional 
energy loss significantly influences the quenching power 
quantified as the overall magnitude of neutral pion Raa 
at RHIC, but that the shape as a function of pr docs 
not show a strong sensitivity. We emphasize that the de- 
scription of Raa is not enough to prove the consistency 
of a specific energy loss mechanism with data, if assump- 
tions about the medium evolution can be rather freely 
adjusted. Therefore folding the jet energy loss mecha- 
nism with a dynamical evolution model which has been 
well tested to describe soft observables is necessary. 
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